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ABSTRACT Introduction of potential disease vectors into a new geographic area poses health risks 
to local human, livestock, and wildlife populations. It is therefore important to gain understanding 
of the dynamics of these invasions, in particular its sources, modes of spread after the introduction, 
and vectorial potential. We studied the population genetics of Aedes (Fmlaya) japonicus Japonicus 
(Theobald), an Asian mosquito that was recognized for the first time in the United States in 1998. 
We examined patterns of genetic diversity using random amplified polymorphic DN A and sequences 
of ND4 of mtDNA by comparing samples from populations spanning the range of this mosquito in 
Japan (six samples) and the United States (nine samples) as well as specimens intercepted in New 
Zealand in 1999. We found geographically differentiated populations in Japan, indicating limited 
gene flow even on small spatial scales. In the United States, we found evidence of significant genetic 
differentiation between samples from New York, Connecticut, and New Jersey and those from 
mid-Pennsylvania and Maryland. We were unable to pinpoint the source location (s) in Japan, 
although some of the U5. samples are genetically close to samples from south Honshu and western 
Kyushu. Further studies should include samples from Korean populations. Distinct genetic signatures 
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if and where they meet. 
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Invading nonindicehous species cause public health 
problems and major environmental damage; and as 
trade and travel increase the number of unwanted 
introductions is expected to rise (OTA 1993). From a 
medical perspective, some of the most worrisome in¬ 
troductions are mosquitoes or other potential disease 
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vectors. The impact of diseases on formerly naive 
human and wildlife populations due to introduced 
mosquitoes is well-documented (Morgan 1981, Bryan 
1999). Furthermore, newly established associations 
between vectors and local parasites may result in 
deadly combinations (Cheng et aL 1999). There is, 
therefore, the potential for major disea-se epidemics 
fueled by newly established associations of parasites, 
hosts, and vectors. To develop the tools to prevent and 
control such epidemics there is a need to understand 
invasion dynamics. That is, what happens to a popu¬ 
lation shortly after it has become isolated from other 
populations of the same species in an unfamiliar and 
possibly ecologically distinct area. 

In February 1999, EL Peyton identified four mos¬ 
quitoes collected in Suffolk County, NY, and Ocean 
County, NJ, in late summer of1998, as Aedes (Finlaya) 
japonicus japonicus (Theobald), establishing a first 
record for the species in the United States (Peyton ct 
al. 1999). Subsequently, Munstermann and Andreadis 
(1999) reported Ae. j. japonicus in Connecticut col¬ 
lected in July 1998. In the spring and summer of 1999, 
large munbers of larvae or adults of Ae. j. japonicus 
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were collected in New Jersey, New York, and Con¬ 
necticut, and a few eggs and larvae were collected in 
Pennsylvania and Ohio. In June of 2000,, Ae. japonicus 
was also found in Frederick, MD (M.R.S., unpublished 
data). 

Aedes j. japonicus is considered one of four mor¬ 
phologically similar subspecies that occur throughout 
most of Japan, Taiwan, Korea, eastern China, and 
Russia (Tanaka et al. 1979). Ae.j. japonicus is relatively 
common in Palearctic Japan (Hokkaido, Honshu, 
Shikoku, Kyushu, Yakushima, and Tsushima) and 
Korea (Peninsula and Cheju Do). Ae. j. japonicus 
is known to feed on chickens and mice (Miyagi 1972) 
and although it is usually found in forested areas it 
win readily bite if humans encroach on its habitat 
(Knight 1969). Ae. j. japonicus breeds in a variety of 
natural and artificial containers very much like the 
better-known Ae. (Stegomyia) albopictus (Skuse). 
This has lead to the suspicion ths-t Ae.j. japonicus, like 
the above-mentioned species, was introduced into the 
United States by the used tire trade (Peyton et al. 
1999). 

Unlike Ae. aegypti and to some degree Ae. albopic¬ 
tus, which are found predominantly in warm climates, 
Ae. j. japonicus is a species adapted to colder condi- 
tions and capable of surviving snowy winters (Tanaka 
et aL 1979). Although Ae.j. japonicus is not considered 
an important disease vector in Japan and Korea, 
laboratory experiments have shown that it is capable 
of transmitting Japanese encephalitis to mice 
(Takashima and Rosen 1959). Because of differeiices 
between habitats, changes in predators, and possible 
alterations in the genetic makeup of the populations as 
a result of the introduction, the potential for Ae. j. 
japonicus becoming a major pest and disease vector in 
the northeastern forested suburbs or even urban areas 
should not be underestimated Recent studies have 
demonstrated that Ae. j. japonicus is a good laboratory 
vector of West Nile virus (M. J. Turell, U.S. Army 
Medical Research Institute of Infectious Diseases, 
Fort Detrick, MD, personal communication), an en¬ 
cephalitis recognized for the first time in the United 
States in 1999 and responsible for several deaths in 
New York City (Sarnpsoii et al. 2000). Indeed, pools ot 
field-collected Ae. j. japonicus have been found posi¬ 
tive for West Nile virus in New York (http://www. 
health.state.ny.us/nysdoh/westnile/index.htm). 

Our objectives were to infer the type of introduc¬ 
tion (unique or multiple) and the geographical origin 
of such introduction (s). Information such as this will 
aid in the development of effective means of control¬ 
ling further introductions. We also aimed to examine 
patterns of genetic diversity among the currently 
known U.S. populations so we can begin to develop 
predictive models of the expansion of this, and possi¬ 
bly other, nonindigenous species. Ultimately, our ob¬ 
jectives are to predict patterns of spread of nonindig¬ 
enous species and tlieir potential roles in the 
occurrence and spread of emerging diseases. 


Materials and Methods 

Included in this study are specimens of Ae. j. ja¬ 
ponicus from ail U.S. populations sampled in 1999 and 
a population sampled in 2000 in Frederick, MD (Table 
1). We also include samples from populations span¬ 
ning the full range of the distribution of Ae.j. japonicus 
in Japan (Table 1), as well as Ae. j. japonicus found in 
a bucket of water on a cargo ship docked in Auckland, 
New Zealand. This species has been intercepted sev¬ 
eral times in New Zealand (Laird et al. 1994) but is not 
established there. Finally, as an out-group we used 
specimens of Ae. j. yaeyamensis from the Iriomote 
Islands (Table 1). This subspecies only occurs in the 
Yaeyama Retto in the Ryukyu Archipelago. Recent 
work on the Ae. japonicus species complex shows tliat 
of the three other subspecies Ae. j. yaeyamensis is the 
closest genetically to Ae. j. japonicus but has a unique 
genetic signature (D.M.F., unpublished data), making 
it suitable as an out-group. We deposited voucher 
specimens in the National Museum of Natural History, 
Smithsonian Institution (accession #1701,1703-1705, 
1708-1709). 

Because all populations of an introduced species 
may be recently derived from the same genetic stock, 
to detect unique genetic signatures among popula¬ 
tions it is necessary to use genetic markers that have 
either very high rates of mutation or were highly 
polymorphic in the original population. Therefore, as 
a first step we chose to examine the genetic diversity 
of Ae.j. japonicus using random amplified polymorphic 
DNA (R.APD). This method js relatively easy and fast 
to implement (Kambhampati et aL 1992), and it allows 
the simultaneous probing of the entire genome. RAPD 
bands can be highly polymorphic because the tech¬ 
nique screens coding and noncoding areas of the ge¬ 
nome (Williams et al. 1990), the latter of which may 
include middle or highly repetitive DNA sequences 
with high mutation rates (Williams et al. 1990, Kazan 
et aL 1993). Unfortunately, analysis of RAPD fre¬ 
quency distributions may also underestimate genetic 
distances and predicted rates of gene flow between 
populations if mutations frequently result in the in¬ 
dependent loss or gain of RAPD bands in different 
populations (Gorrochotegui-Escalante et al. 2000). 
This problem should, however, be less important 
when examining recently established populations like 
those resulting from an introduction. 

As a second class of marker, we sequenced 424 bp 
of ND4, a mitochondrial gene (NADH dehydrogenase 
subunit 4). mtDNA is maternally inherited (Avise 
1994) and mtDNA variants segregate rapidly between 
generations leading to relatively high rates of poly¬ 
morphism (Howell et al. 1996, Jenuth et aL 1996). 
From the sequence data one can recover both rates of 
gene flow and phylogenetic relationships between 
populations. For these reasons, mtDNA has been used 
extensively in population genetic studies (e.g.. Conn 
et al. 1993, Besansky et al. 1997, Fonseca et al. 2000). 
In particular the ND5 and ND4 mtDNA loci have been 
shown to be highly polymorphic in both Anopheles and 
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Table 1. Atdet japonicus samples, list of locations with latitude and longitude, date and Ufe stage at which specimens wete collected, 
and sample sizes for RAPD ami ND4-mtI>NA analyses 



Aedes species (Lehmann et al. 1997, De Merida et al. 
1999, Gorrochotegui-Escalante et al. 2000). 

Because we only had four to six individuals from 
Ohio, Pennsylvania, Maryland, and New Zealand, we 
chose to examine a minimum of four and a maximiun 
of six specimens from all locations and we favored 
analyses that treat specimens individually. Once we 
established the degree of similarity between individ¬ 
uals we combined specimens from several locations to 
achieve sample sizes large enough to support popu¬ 
lation level analyses (Nei 1978). When we had large 
samples (New York, New Jersey, Connecticut, and 
Japanese locations), we chose specimens from as 
many different soiurces and collecting trips within a 
location as possible (e.g., different dates, water cavi¬ 
ties, or pools of mosquitoes). The specimens collected 
in Chitose, a city in Hokkaido, Japan, and those in¬ 
tercepted in Auckland, New Zealand, did not amplify 
well with the RAPD primers possibly because of DNA 
degradation during transport to the United States. 
Specimens from several locations in Connecticut 
where improperly preserved and had developed fungi. 
For that reason we excluded them from the RAPD 
analysis. Since we were able to amplify ND4 success¬ 
fully from many of these specimens, however, they 
were included in the ND4 analysis. Because there 
were only one to three specimens from Essex, Port¬ 
land, Roxbury, or Stamford we combined them with 
the foiu specimens from Kent to form a "Connecticut” 
sample (Table 1). 

RAPD Loci. We extracted DNA from individual 
specimens using a phenol/chloroform method (Fon¬ 
seca et al. 2000). We screened 30 decamer primers 
(AOl-07,09,10,12-16,20; BOl-05.16-18; C03,04,13, 


18,19; Dll, 13, Operon Technologies, Alameda, CA), 
and used nine in the analyses. Procedures followed 
very closely those described in Wilkerson et aL (1995) 
except that we used a PE Biosystems 9700 thermocy- 
cier (PE Biosystems, Foster City, CA) and we ran the 
gels for 5 h. The polymerase chain reaction (PCR) 
parameters were 5-min denaturation at 94°C followed 
by 45 cycles of 1 min at 94°C, 1 min at 35°C, and 2 min 
at 72°C. Minimum ramp times were used at all stages 
except between 35 and 72°C, where we delayed the 
ramp by 25% to achieve a rate of 4°C/s. Molecidar 
weight standards to size the RAPD bands were mostly 
provided by phiX174 DNA digested with HaeUI 
(Sigma, St. Louis, MO). In the few instances in which 
the RAPD bands were larger than 1,350 bp (primers 
A20 and C13), we used a 3;1 mix of phiXI 74Hae m and 
lambda DNA digested with ffindin (Sigma). The mo¬ 
lecular weight of each RAPD band was estimated by 
comparison with the co-migrating size standards using 
a least-squares polynomial regression analysis. For 
each primer a minimum of 15 randomly chosen sam¬ 
ples were rerun at least once, in some cases two to 
three times. Only bands shown to be reproducible, and 
in the intermediate molecular weight range, were 
scored Because the preservation of specimens dif¬ 
fered somewhat across samples, the choice of primers 
was based on their ability to consistently produce 
bands of equivalent size range and intensity across all 
the samples examined. Although this may result in an 
underestimate of differences between samples, we felt 
the known sensitivity of the RAPD method to DNA 
quality (Koch et al. 1998) made such a concession 
necessary. We also examined specimens carefully un¬ 
der a dissecting microscope before DNA extraction 
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and discarded the entire sample for RAPD analysis if 
any specimens had fungi or signs of possible contam¬ 
ination with other sources of DNA. 

The RAPD bands were analyzed as genetic markers 
assuming that recessive “band absent” alleles are iden¬ 
tical among and within individuals, and dominant 
“band present” alleles are identical among and within 
individuals (Gorrochotegui-Escalante et al. 2000). 
When performing exact tests we also assumed that 
RAPD bands of different sizes result from indepen¬ 
dent mutations. The analyses were performed either 
by considering specimens individually or by combin¬ 
ing them into groups. For the analysis of individuals 
the binary data matrix was converted to a similarity 
matrix using Nei and Li (1985) similarity index and 
then to a distance matrix (D = 1-S), in RAPDPLOT. 
The statistical package RAPDPLOT was developed by 
William C. Black and is available by anonymous ftp 
(ftp://lamar.colostate.edu/pub/wcb4/). Specimens 
were combined into groups based on the individual 
analysis. We usedTFPGA version 1.3 (Miller 1997) to 
calculate Nei’s (1978) unbiased heterozygosity and 
the percent of polymorphic RAPD loci in each sample. 
We also used TFPGA to calculate exact tests of pop¬ 
ulation differentiation (Raymond and Rousset 1995). 
RAPDDIST in RAPDPLOT was used to compute pair¬ 
wise distances between samples. We analyzed the data 
using all the genetic distance measures available in 
RAPDDIST with the correction Lynch and Milligan 
(1994) proposed for dominant markers like those de¬ 
rived from RAPD analyses. As a graphic representa¬ 
tion of the relationships among specimens and sam¬ 
ples, dendrograms were constructed by unweighted 
pair group method with arithmetic averaging (un¬ 
weighted pair-group method with arithmetic average, 
UPGMA Swofford et aL 1996) using the NEIGHBOR 
and GONSENSUS programs from PHYLIP version 
3.57c (Felsenstein 1995). To evaluate the relative 
strength of the tree nodes we analyzed 1,000 bootstrap 
replicate matrices generated with RAPDPLOT. Trees 
were visuahzed with TreeViewPPG (Page 1996). We 
tested for the possibility of isolation by distance (Slat- 
Idn 1993) by comparing distance matrices with a ma¬ 
trix of geographic distances between locations in Ja¬ 
pan and the United States with a Mantel test (Mantel 
and Valand 1970) in Arlequin version 2000 (10,000 
permutations, Schneider et al. 2000). 

ND4 Sequencing. Based on the universal primers 
reported by Simon et al (1994) and the complete mi¬ 
tochondrial sequence for Anopheles gambiae (Beard et 
al. 1993) and An. cpiadrimaculalm (Mitchell et aL 1993) 
we designed two new primers; N4J-8502D 5'-GGTAG- 
GAGGAGCAGCTATATT-3' and N4N-8944D 5'-AAG- 
GCTCATGTTGAAGCrrGG-3'. These primers amplify 
a 424-bp fragment between positions 8398 and 8821 in 
the An. gambiae sequence (GenBank accession 
#L20934). For the amplification we used 5 ng of 
genomic DNA in a 50-p.l reaction. The final concen¬ 
trations of the PGR reagents were as follows; 1 X PGR 
Buffer, 300 nM of each primer, 250 (j,M of each dNTP, 
2 mM MgC12, and 1.5 U of Tag Gold polymerase (PE 
Biosystems). The PGR amplification was preceded by 


a 10-min denaturation at 96°G. The amplification con¬ 
sisted of 35 cycles of 40 s at 94°G, 40 s at 55°G, and 60 s 
at 72°C, and ended with a final extension step of 7 min 
at 72‘’G. Some individuals that did not amplify under 
these conditions were subsequently amphfied by low¬ 
ering the annealing temperature to SS'C. After clean¬ 
ing the PGR products witli QiaQuick columns (Qia- 
gen, Valencia, GA), cycle sequencing was done in both 
directions using the PGR primers and ABI Prism Big- 
Dye Terminators (PE Biosystems). We used one-half 
of the 1:4 dilution conditions described in PE Biosys¬ 
tems Application Note #107PAB01-01 in a total vol¬ 
ume of 10 pi. Sequences were visuahzed and scored on 
a 373 PE Biosystems Automatic Sequencer (PE Bio¬ 
systems). 

Alignment and assembly of sequences were per¬ 
formed with Sequencher 3.0 (Gene Godes, Ann Har¬ 
bor, MI). Arlequin version 2.000 (Schneider et al. 
2000) was used to calculate haplotype frequencies, 
gene diversity (H) and its samphng variance, and 
mean number and variance of pairwise differences 
between haplotypes in each sample (■jr). H is equiv¬ 
alent to the expected heterozygosity for diploid data 
(Nei 1987) and is defined as the probabifity that two 
randomly chosen haplotypes are different in a sample. 
The total variance of rr is an estimate over the sto¬ 
chastic and the sampling processes (Tajima 1983). 
After these initial analyses by sample, we combined 
the specimens into groups based on putative haplo¬ 
type discontinuities and used Arlequin to estimate 
pairwise values from haplotype frequencies and 
perform exact tests of population differentiation. Be¬ 
cause of the large number of tests performed we used 
a sequential Bonfeironi correction to adjust the a 
value for each test (Rice 1989). Basic statistical tests 
were performed with JMP version 3.1.6 (SAS Institute 
1995). 

Results 

RAPD Analyses. The nine primers used are AOl-03, 
A20, B17, B18, G13, G18, and Dll and together they 
yielded a total of 112 RAPD bands (93 bands in anal¬ 
yses that did not require an out-group) after we re¬ 
moved from the analyses those bands that occurred in 
only one individuaL Only three bands were larger than 
1,350 bp (from primers A20 and C13) and only two 
were smaller than 310 bp (both from Dll, found ex¬ 
clusively in the Pennsylvania and Maryland samples). 
The presence/absence matrix is available (http:// 
wrbu.si.edu/wrbu.html). The trees that resulted from 
the UPGMA analysis did not vary with the type of 
genetic distance being used although occasionally a 
bootstrap value changed from <50% to >50%, de¬ 
pending on the analysis. In these cases we report the 
lower value. There were consistently only three 
RAPD bands in common between Ae. j. japonicus and 
Ae. j. yaeyamensis, which underlines the genetic dis¬ 
tinctness of the latter subspecies and its usefulness as 
an out-group in these analyses. 

Analysis of the five Japanese samples revealed sig¬ 
nificant differences between locations. All specimens 
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from Sapporo, Tokyo, Nagasaki, and most of those 
from Saga and Hiroshima formed separate groups 
(Fig. 1). One specimen from Saga and one from 
Hiroshima, however, did not cluster within their re¬ 
spective locations (Fig. 1). Either considering Saga 
and Hiroshima as separate populations or not, pairwise 
measures of genetic distance are significantly corre¬ 
lated to geographic distance (r^ = 0.40, Mantel prob- 
ahility = 0.041 for Nei’s (1978) unbiased genetic dis¬ 
tance with Saga and Hiroshima as separate samples). 

In the United States, we found an association be¬ 
tween the specimens collected in New York, Con¬ 
necticut, and northern New Jersey (yemon) loca¬ 
tions (Fig. 2). The specimens from New Egypt, NJ, 
form a separate cluster, as do those from Oak Hill, OH, 
but Pennsylvania and Maryland samples cluster to¬ 
gether (Fig. 2). Still, the degree of differentiation 
between locations is considerably less than that found 
in the analysis of the Japanese samples (Fig. 1) and we 
found no correlation between genetic distance and 
geographic distance (r* = 0.02, Mantel probability = 
0.32). 

ND4 Analysis. We sequenced this locus for 82 spec¬ 
imens of Ae. j. japonicus, including those used in the 


RAPD analysis. We found 15 haplotypes that result 
from 15 transitions (Fig. 3, CenBank accession 
#AF305879). Five haplotypes were only found in 
Japan, one is unique to New Zealand, and six are 
unique to the United States (Fig. 3). Only one hap- 
lotype, HI, is common to the United States, Japan, and 
New Zealand. On average we found 2.3 ±0.3 (mean ± 
1 SE) ND4 haplotypes in each location in Japan and 
2.2 ± 0.2 in the United States analysis of variance 
(ANOVA; P = 0.78). Three haplotypes were found in 
the four specimens obtained in Auckland, New Zea¬ 
land (Fig. 4). The gene diversity is not significantly 
different between Japanese and U.S. samples (0.61 ± 
0.10 in Japan and 0.47 ± 0.85 in the United States, 
ANOVA, F = 0.99; df = 8,5; P = 0.20) nor is the mean 
munber of pairwise differences between haplotypes 
(1.53 ± 0.36 in Japan and 1.38 ± 0.29 in the U.S. 
samples, ANOVA, F = 0.11; df = 8, 5; P = 0.44). 

In Japan, the patterns of ND4 diversity mirror those 
found with RAPDs. Samples from Hokkaido, the 
northernmost island in Japan (Sapporo and Chitose), 
have unique haplotypes although the ubiquitous hap- 
lotype HI was found in the Chitose sample. Likewise, 
Tokyo specimens have unique haplotypes, although 
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the presence of haplotype H3 indicates possible ge¬ 
netic exchange with the southwestern populations. In 
the United States, two patterns are evident. First, tlie 
five samples of populations clustered cn and around 
Long Island (Shelter Island, Southold, Manorville, 
New Egypt, and Vemon, in New Jersey and the five 
locations in Connecticut) all appear to be successive 
samplings from the same gene pool. The highest ge¬ 
netic diversity occurs in Long Island, where five hap- 
lotypes can be found (Table 2). Second, Pennsylvania 
(PA) and Maryland (MD) have one unique haplo¬ 
type, H9, which occurs with high frequency (Fig. 4) 
and a second haplotype, H12, also found nowhere else 
in the United States, although it occurs in one indi¬ 
vidual from Chitose in Northern Japan (Fig 4). Ohio 
has one unique haplotype, H4, but the most common 
haplotype there is HI. If present, haplotypes H4, H9, 
and Hli must occur with frequencies lower than 0.1 
in the New York, Connecticut, and New Jersey pop¬ 
ulations because, after sampling 35 individuals, they 


have stiU not been detected (the probability of not 
samphng a haplotype if that haplotype is present = 
[ 1-frequency] where N is the number of specimens 
tested). 

Combined Analysis. Based on the analyses done 
separately for Japanese and US. samples, we com¬ 
bined specimens into eight groups (Table 2) and per¬ 
formed a population level analysis with all specimens 
and the out-group. The RAPD analysis shows a pattern 
of interspersion (Fig. 5). The NY/Veraon/CT group 
is more closely related to specimens from southwest¬ 
ern Japan (Hiroshima and Saga) than to other U.S. 
groups. Significant Fst values and exact tests of pop¬ 
ulation differentiation with the ND4 and RAPD data, 
respectively (Table 3), show differentiation between 
NY/Vemon/CT and Pennsylvania/Maryland as well 
asbetween NY/Veraon/CT and Hokkaido samples. At 
the mtDNA-ND4 locus, 68% of the pairs of samples 
analyzed were stcitistically different. 
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Fig. 3. Mitochondrial (ND4) DNA haplotypes found in this study and their observed frequencies. Only variable loci are 
shown and the numbers above each represent its position along the 434-bp length of the ND4 gene sequenced. Haplotypes 
found in the United States (except for Pennsylvania and Maryland) are shown in black, those found in Japan are shown in 
white. Haplotypes from the Pennsylvania/Maryland samples are in gray and haplotypes from New Zealand are in bght 
stippled. 


Discussion 

In Japan, both RAPD and mtDNA analyses reveal 
differentiation among populations. We found unique 
RAPD banding patterns and mtDNA hsqilotypes in 
most locations and a systematic auid significant de¬ 
crease in genetic similarity with distance between 
locations, a signatinre of long-established populations 
isolated by distance (Rousset 1997). Because Ae. ja¬ 
ponicus specimens have heen able to arrive in the 
United States and were found in cargo ships in New 
Zealand, we were surprised by not finding extensive 
gene flow between Japanese locations. Although 
movement across Japan might be precluded by its 
insular nature, there were significant differences even 
between locations in the same island (e.g., Hiroshima 
and Tokyo). The distribution of Ae. japonicus in Japan 
suggests therefore that this mosquito has a relatively 
short flight range and is not commonly transported by 
humans. Still, the presence of the mtDNA haplotype 
HI in the saunple from Chitose in northern Japan could 
indicate that some human mediated gene flow be¬ 
tween Japanese locations is occurring. Based on a 
minimum spanning tree analysis (data not shown), 
however, HI is putatively the haplotype from which 
50% of the remaining haplotypes evolved. It is then 
possible that the high frequency of HI across most 
populations reflects not contemporaneous gene flow 
but an ancestral common source of all populations. 
iTiat would also explain its abundance in the U.S. 
samples, although most of the remaining haplotypes 
are unique to either the U.S. or Japanese populations. 
It should be noted that the presence of haplotype H5 
in the New Zealand sample, a haplotype only found in 
Hokkaido samples, does appear to indicate a northern 
Japanese source for those mosquitoes. Therefore, 


some human-associated traffic of Ae. japonicus from 
Japan is occurring. 

In the United States, different RAPD patterns and 
mtDNA haplotypes, particularly between NY/NJ/CT 
and Pennsylvania/Maryland specimens, but also the 

x> A _ __ VT —*. kit j 

Oak Hill, OH, might suggest multiple introductions. 
Because the sample sizes in these analyses, the first 
performed on this species, were constrained by the 
small number of specimens available, the power (Co¬ 
hen 1988) of the analysis is relatively low. Neverthe¬ 
less, it is undeniable that there are indications of at 
least two independent foci of expansion of Ae. j. ja¬ 
ponicus in the USA, one located in the northeastern 
states of New York, Connecticut, and New Jersey, and 
the second in the main axis of Pennsylvania and Mary¬ 
land. These two foci are significantly genetically dif¬ 
ferentiated as evidenced by the significant Fst values 
and exact tests for population differentiation (Table 
3). Although ND4-mtDNA haplotypes are very dis¬ 
tinct between New York/New Jersey/Connecticut 
and Pennsylvania/Maryland, RAPD patterns between 
these two groups of samples do not appear as distinct 
as those found between locations in Japan. The RAPD 
patterns may therefore reflect a common source for all 
U.S. populations that, once isolated, will develop very 
different ND4-mtDNA haplotype frequencies. On the 
one hand, the fact that a single gravid mosquito may 
start a new population makes such a scenario feasible. 
On the other hand, the virtual absence of haplotypes 
H9 and H12 from the New York/New Jersey/Con- 
necticut samples, makes the coincidence of two gravid 
mosquitoes (each female carries a single mtDNA hap¬ 
lotype) , with very rare haplotypes founding the Penn¬ 
sylvania/Maryland population, an unlikely scenario at 
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Fig. 4. MtDNA-ND4 haplotype distribution and geographic location of all populations sampled for this study, (a) U.S. 
locations, (b) Japanese locations. Numbers in the pie charts correspond to column one of Table 1. Each color represents a 
different haplotype (H1-H15, Fig. 3), and in each pie chart the area of a color corresponds to the proportion of the respective 
haplotyiie in that sample. 


best. It remains possible, however, that both New 
York/New Jersey/Connecticut and Pennsylvania/ 
Maryland populations derived from an original intro¬ 
duction in which all U.S. ND4-mtDNA haplotypes 
were relatively common. Such a population would 
have an mtDNA genetic diversity similar to that en¬ 
countered in some Japanese locations (e.g., Tokyo or 
Chitose) and it still might be present in the United 
States, undetected or at least untested. 

After a founder event like a new introduction, one 
might expect to find reduced genetic diversity (Haiti 


and Clark 1997). Instead, neither RAPD nor ND4- 
intDNA analyses show evidence of an overall reduced 
genetic diversity in the United States when compared 
with Japanese samples (Table 2). Although the sam¬ 
ples from Pennsylvania and Maryland have a signifi¬ 
cantly lower ND4-mtDNA diversity than samples 
from southern Japan (ANOVA, F = 12.57; df = 9,19; 
P = 0.02) this pattern is not mirrored in the RAPD 
band diversity. There are at least three possible ex¬ 
planations for the high genetic diversity found in U.S. 
populations ofAe.japonicus: (1) the species arrived in 
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Table 2* Measures of genetic variation in population level analyses using both RAPDs and mtDNA-ND4 




RAPDs 




ND4 


Locations 

n 

H 

%p 

n 

#h 

H £ V(H) 

•jr± V(7r) 

Southold, NY 

5 

0.12 

33.3 

5 

2 

0.60 ± 0.18 

1.20 ± 1.07 

Shelter Island, NY 

5 

0.13 

34-4 

5 

2 

0.40 i 0.24 

1.20 ± 240 

Manorville, NY 

5 

0.12 

32.1 

5 

3 

0.70 ± 0.22 

2.20 ± 307 

Connecticut 

4 

0.11 

29.0 

10 

3 

0.51 ± 0.16 

0.91 ± 1.00 

Vernon, NJ 

All Long Island locations, 

5 

0.11 

26.9 

5 

3 

0.70 ± 0.22 

1.20 ± 1.07 

Vernon (NJ), Connecticut** 

24 

0.15 

55.9 

30 

5 

0.73 ± 0.05 

1.55 ± 1.42 

Ohambersbui^, PA 

5 

0.13 

34.4 

6 

2 

0.33 ± 0.22 

1.67 ± 5.95 

Frederick, MD 

5 

0.11 

30.1 

4 

2 

0.50 ± 0.27 

2.50 ± 7.50 

Pennsylvania+Maryland* 

10 

0.14 

43.0 

10 

2 

0.36 ±0.16 

1.78 ± 5.86 

New Egyjrf, NJ“ 

5 

0.09 

22.6 

5 

1 

0.00 i 0.00 

0.00 ± 0.00 

Oak Hill, OH' 

4 

0.11 

25.8 

4 

2 

0.50 ± 0.27 

1.50 ± 2.70 

All US locations 

43 

0.17 

75.3 

49 

8 

0.83 ± 0.03 

2.41 ± 2.55 

Nagasaki* 

4 

0-10 

24.7 

5 

2 

0.60 ± 0.18 

240 ± 4.27 

Saga 

4 

0.14 

37.6 

5 

2 

0.70 ± 0.22 

200 ± 4.00 

Hiroshima 

4 

0.13 

33.3 

5 

3 

0.70 ± 0.22 

0.80 ± 0.40 

Saga+Hiroshima' 

8 

0.15 

51.6 

10 

4 

0.73 ± 0.10 

1.56 ± 2.89 

Tokyo* 

4 

0.16 

37.6 

5 

3 

0-80 ± 0.16 

2.80 ± 2.40 

Sapporo* 

4 

0.12 

30.1 

5 

1 

0.00 ± 0.00 

0.00 ± 0.00 

Chitose* 

— 

— 

— 

4 

3 

0.83 ± 0.22 

1.17 ± 0.57 

All Japanese locations 

20 

0.19 

77.4 

29 

8 

0-84 ± 0.04 

2.13 ± 2.71 

Auckland, NZ 

— 

— 

— 

4 

3 

0.83 ± 0.22 

1.67 ± 1.47 


n, sample size for each analysis; H, Nei’s (1987) unbiased heterozygosity; %P, percent RAPD loci not fixed; #h, number of haplotypes; H, 
gene diveisity; V(H), sampling variance; ir, mean number of pairwise differences between haplotypes; V(7r), total variance in mean number 
of pairwise differences. 

' Samples were used to examine the similarity between US and Japanese populations. (Fig. 5). 


large numbers, (2) there were multiple introductions, 
and (3) Ae. japonicus has been in the United States for 
quite some time, undetected. Black et aL (1988), 
Kambhampati et al. (1990) andUrbanelJi et al. (2000) 
argued that large introductions might explain why 
they failed to find reduced genetic diversity in U.S. and 


Italian populations of Ae. albopictus. The two most 
common scenarios for the arrival of mosquitoes are as 
adults in passenger planes or as desiccated eggs in tires 
or other similar containers in cargo ships. The number 
of mosquitoes that arrive on a plane is almost certainly 
small (Raymond et al. 1991, Qiao and Raymond 1995) 



H-1--1 

0 ^ 0.1 0.0 

Genetic distance 


Fig. 5. UPGMA dendrogram based on RAPD loci using Nei’s standard genetic distance showing putative relationships 
among groups of Ae, j. japonicus specimens from all the locations included in this study. fXitgroup is Ae.J. yaeyamensis. Sample 
sizes (N) are shown in front of each group label. Bootstrap support is based on 1,000 permutations. Only bootstrap values 
higher than 50% are shown. 
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TaUe 3. Pairwise Fst ▼ahies between populations and exact tests for population differentiation 



1 

2 

3 

4 

5 

6 

7 

8 

1 


0.301 

0,303 

0.000" 

0.636 

0.375 

0.038 

0.000" 

2 

0.545“ 


0.794 

0.140 

0.984 

0.992 

0.753 

0.192 

3 

0.092 

0.761" 


0.834 

1.000 

1.000 

1.000 

1.000 

4 

0.552" 

0.756" 

0.389" 


0.924 

0.201 

0.603 

0.015 

5 

0.199 

0.665" 

-0.013 

0.498" 


1.000 

1-000 

0.743 

6 

0.278“ 

o.eer 

-0.094 

O.ITT 

0.161 


1.000 

0.984 

7 

0.258" 

0.533 

0.027 

0.452“ 

0.156 

0.174" 


1.000 

8 

0.179* 

0R62" 

0.433 

0.634" 

0.470" 

0.419" 

0-448“ 



Values below and above the midline are Fs^- values from the mtDNA-ND4 analysis and P-values from exact tests for population differentiation 
from the RAPD data, respectively. 1, NY/ Vemon/ CT; 2, New Egypt; 3, Ohio; 4, PA/ MD; 5. Nagasaki; 6, Saga+Hiroshima; 7, Tokyo; 8, Sapporo 
(RAPD) or Sapporo+Chitose {ND4). Sample sizes as in Table 2. 

* Significant value using sequential fionferroni corrected ct values. 


but may happen repeatedly, whereas eggs in tires 
could result in large populations if the tires get 
flooded, although, the effective population size may 
be considerably smaller. Of notice is the fact that 
although the tire trade has been implicated in the 
introduction of many species of mosquitoes to many 
places (Craven et al. 1988, Reiter 1998), extensive 
examinations of tire transport into the United States 
did not reveal the presence of Ae. japonicus. Laird et 
al. (1994) found Ac. japonicus in only three out of 8,549 
tires inspected that arrived in New Zealand between 
November 1992 and January 1993. Overall, therefore, 
Ae. japonicus might not be as common an occurrence 
in tires as Ae. albopictus. 

How Ae. japonicus is moving between locations 
within the United States is also not clear. Certainly the 
tire trade within the United States might be an im¬ 
portant source of new focL Another possibility is that 
mosquito movement is associated with the Standard- 
bred horse trade. Suffolk County, NY, and New Egypt, 
N J, where Ae. japonicus were first collected, are breed¬ 
ing centers of Standardbred horses. The location in 
Vemon, NJ, where Ae. japonicus were collected in 
large numbers is also near a Standardbred farm. Adult 
Ae. japonicus may use Standardbread trailers as resting 
areas and be distributed between racing areas and 
farms. If rare, this form of distribution might decrease 
the vigor of the Ac. japonicus colonization of the 
United States and increase the effectiveness of control 
measures. Indeed, it does appear that Ae. japonicus is 
cnrea din ff more slowly than Ae. albopictus did. In 1987, 
“year two” after its detection in large numbers in 
Texas, Ae. albopictus was already known from 12 sep¬ 
arate states (Moore et al. 1988). In contrast, by July 
2000 thenumber of states with confirmed Ae. japonicus 
records has only increased from two to six. It should 
be noted, however, that the first U.S. Ae. albopictus was 
actually caught in Tennessee in 1983. That introduc¬ 
tion was traced to the used-tire trade but thought to 
be an isolated incident (Reiter and Darsie 1984). Fur¬ 
thermore, Ae. alhopictus had been identified in tires 
arriving in American ports as early as 1946 (Pratt et al. 
1946, Eads 1972). It is then possible that by 1987, Ac. 
albopictus had already been multiplying in the United 
States for several years, undetected. We suspect that 
the extensive network of Vector Control Agencies and 
protocols established after the introduction of Ae. al¬ 


bopictus to the continental United States (Moore 1986, 
Moore et al. 1988, Francy et al. 1990) allowed the 
identification of Ac. japonicus soon after its arrival. 
Still, the evidence of high levels of genetic diversity of 
U.S. Ae. japonicus populations suggests that better 
measures have to be implemented to control the im¬ 
portation of exotic species. 

We did not find a perfect match in Japan for any of 
the U.S. samples. Ejccepting the putatively ancestral 
and very common ND4-mtDNA HI haplotype and the 
presence of haplotype H12 in the sample from Chi- 
tose, which warrants a broader sampling effort in 
northern Japan, no other haplotype was found in com¬ 
mon between U.S. and Japanese samples. Although 
the high diversity of mtDNA-ND4 haplotypes in Japan 
suggests we may simply have missed sampling the 
source populations, the RAPD band patterns correlate 
the NY/NJ/CT samples closest with southwestern 
Japanese locations, which we sampled repeatedly. Be¬ 
cause these areas have extensive commercial contact 
with South Korea (D. Strickman, Walter Reed Army 
Institute of Research, Silver Spring, MD, personal 
communication), it is clear that we need to do a 
broader reconnaissance of the genetic diversity of this 
species, including analysis of Korean populations, to 
address the question of source and means of trans¬ 
portation. 

Our study shows at least two genetically indepen¬ 
dent foci of expansion of Ae. japonicus in the United 
States. It is unclear at this point if they originated from 
separate introductions or if populations have diverged 
since they have arrived. More extensive sampling and 
testing of U.S. populations will be needed to address 
this question. If several introductions are involved, 
their present location has placed distinct genotypes 
close enough that hybridization is now most probable. 
Hybridization might result in localized increases in 
genetic diversity and in the generation of new genetic 
combinations with unknown vectorial capacity. Al¬ 
though when compared with Ac. j'. yaeyamensis, Ae.j. 
japonicus appears to be a single and distinct genetic 
unit, we have identified discontinuities that may be 
correlated to variation in vectorial potentiaL Because 
we may have been able to examine the genetic diver¬ 
sity of Ae. japonicus early in its introduction history to 
the United States we may also be able to pinpoint 
populations to target for control We will be pursuing 
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our work on the genetics of this group both to under¬ 
stand its impiortance as a vector of human disease and 
as a model of mosquito invasion dynamics. 
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